YITP-11-61 



Large and strong scale dependent bispectrum in single field inflation from a sharp 

feature in the mass 

Frederico Arroja 1 *, Antonio Enea Romano 2 ^ and Misao Sasaki 3 -' 1 
1 Institute for the Early Universe, Ewha Womans University, Seoul 120-750, Republic of Korea 
2 Leung Center for Cosmology and Particle Astrophysics, 
National Taiwan University, Taipei 106, Taiwan 
3 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, 
Japan; Korea Institute for Advanced Study 207-43 Cheongnyangni 2-dong, 
Dongdaemun-gu, Seoul 130-722, Republic of Korea 
(Dated: July 14, 2011) 

We study an inflationary model driven by a single minimally coupled standard kinetic term scalar 
field with a step in its mass modeled by an Heaviside step function. We present an analytical 
approximation for the mode function of the curvature perturbation, obtain the power spectrum 
analytically and compare it with the numerical result. We show that, after the scale set by the step, 
i— ( ' the spectrum contains damped oscillations that are well described by our analytical approximation. 

We also compute the dominant contribution to the bispectrum in the equilateral and the squeezed 
limits and find new shapes. In the equilateral and squeezed limits the bispectrum oscillates and it 
C^) ■ has a linear growth envelope towards smaller scales. The bispectrum size can be large depending 

' on the model parameters. 

I. INTRODUCTION 

u: 

The inflationary paradigm solves successfully the problems of the Big Bang model and also provides us 

Oh' with a mechanism to generate the primordial perturbations from quantum fluctuations that are later imprinted as 
anisotropies in the cosmic microwave background (CMB) radiation and in the large-scale structure of the universe. In 
particular, regarding the CMB power spectrum, we now know that it is nearly scale invariant Q, mostly of adiabatic 
origin and it contains almost all the information about the primordial perturbations (i.e. the perturbations are nearly 
. Gaussian). This is in agreement with the predictions of many models of inflation. 

However, a small amount of non-Gaussianity is still allowed by the CMB data and recently there have been several 
, claims of detection of non-Gaussianity of the primordial perturbations [8T4llj]. If the values of the non-Gaussianity 
parameters are of the order of magnitude claimed then NASA's satellite, Planck fl2|. which is already up in the 
sky taking data, should be able to measure these values at many standard deviations [la ]. These are exciting times 
OO . for theoretical cosmologists that for the past ten years or so have been trying to construct models where large non- 
Gaussianity is produced. 

It is well known that the simplest and most popular inflationary model, a single scalar field with standard kinetic 
term satisfying the slow-roll conditions and with standard initial conditions for the vacuum state of the quantum 
perturbation predicts a level of non-Gaussianity that is small and unobservable fl4l - fl9l ]. even for the Planck satellite 
[20|-[22j. Therefore, the previously mentioned observations have the potential if confirmed to rule out a large and 
most popular class of models. It is good on itself to rule out many models in one go and even better that a detection 
of non-Gaussianity (the same applies to the tensor-to-scalar ratio) would allows us to progress significantly in our 
understanding of the mechanism that drove inflation in the early universe, this is because higher-order statistics 
contain much more information about the dynamics that the power spectrum. 

Our searches for inflationary models producing large observable non-Gaussianity have been fruitful. Many possi- 
bilities have been found, for example models with non-canonical kinetic terms (like DBI-inflation, k-inflation, ghost- 
inflation) [231433 ] , multiple field models of inflation [35l - l66j , temporary violations of the slow- roll conditions and small 
departures of the initial vacuum state from the standard Bunch-Davies vacuum [67T4751 ]. For recent reviews about 
these mechanisms to produce non-Gaussian perturbations see [tgMsoT ] . 

Both at the bispectrum level (three-point correlation function) and the trispectrum level (four-point correlation 
function) many observationally distinct shapes of these higher-order correlations have been found [8ll - [85| . Some 
of these shapes of non-Gaussianity have been constrained with CMB data @, |86| and large-scale structure data 
pll I22I [87-89]. Current limits on the amplitude of the bispectrum, taken from the 7- year WMAP [9(| data at 95% 
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confidence level are Q: -10 < f^ c L al < 74^-214 < f^ a < 266 and -410 < /^' lo < 6 for the local, equilateral 
and orthogonal shapes respectively. Ref. [86| has a more recent analysis of the CMB bispectrum (see [91| for the 
trispectrum analysis) including many other shapes. Planck satellite will improve all these constraints significantly and 
will also improve existing constraints on the amplitude of the trispectrum. For some review papers on observational 
aspects of non-Gaussianity see for instance [92- 96]. 

In this paper, in order to produce large and observable non-Gaussianity, we shall study a model where the slow-roll 
conditions are temporarily violated. This model can be included in the class of mode ls with "features" (i.e. non-scale 
invariant models) in the potential or more generally in the field's Lagrangian [97H103I ]. There are several other reasons 
that motivate the study of "feature(s)" models. For instance, it has been s hown tha t these models can provide better 
fits to the power spectrum of the CMB anisotropics than the ACDM model |104| - |108| . In these models, this is achieved 
thanks to the introduction of new parameters (and scales) that can be tuned to coincide with well known "glitches" 
£ ~ 20 — 40 in the CMB spectrum data . These models are also of theoretical interest because they have very distinct 
bispectrum signatures like scale-dependence (e.g. "ringing" and localization of Jnl)- The features often have a more 
fundamental ph ysica l exp lanation, like for example they might be local ized sharp features due to particle production 
during inflation due to a duality c ascade in brane inflation |ll3j , peri odic features due to the production 

of instantons in axion monodromy inflation jl!4l4ll6l | or phase transitions [l!7j . The study of these more realistic 
scenarios might eventually allows us to identify the underlying microscopic theory of inflation. 

In particular, in this paper we will study a model where the mass of the inflaton field suddenly changes. We will 
describe this using a toy model where the change is approximated by a Heaviside step function. For example, this 
might be a toy model for a first order phase transition. We will compute the power spectrum and the bispectrum of 
the curvature perturbation and show that its amplitude can be large and that its shape is new with very distinctive 
features. 

This paper is divided into the following sections. In section [TTl we introduce the model and give the analytical 
background solution under some approximations. In section Hill we discuss linear perturbations, present an analytical 
approximation for the mode function of the primordial curvature perturbation and finally calculate the analytical 
power spectrum and compare it with the numerical result. In section IPVl we compute the dominant contribution to 
the bispectrum of the curvature perturbation in two interesting limits, namely the equilateral limit and the squeezed 
limit. Section W\ is devoted to the conclusion. 



II. THE MODEL 



We will consider a canonical minimally coupled single scalar field model with the potential given by 

v( ,v _( V Q + \m 2 A 2 : > 0o m 

V(<t>) -\V Qa + \m%l + A)^ :0<0 O (1) 

where (f> is the scalar field, m$ its mass, Vq and Vo a are the vacuum energies before and after the transition field 
value 4>o respectively. A is a constant amplitude and a parameter of the model. We restrict the allowed range of A 
as A > — 1. Requiring continuity of the potential across 4>o implies Vo a = Vo — l/2rr£A<fi 2 l , where <pQ = 0(to) an d 
to is the transition time. We will assume that this change in vacuum energy is small and we shall neglect it. In the 
next section, we will consider perturbations. We will perform the calculation in the comoving gauge therefore the 
background value of <j), i.e. 0o> fully determines when the transition happen even in the perturbed spacetime. 

We are interested in flat, homogeneous, and isotropic Friedmann-Lemaitre- Robertson- Walker background universes 
described by the line element 

ds 2 = -dt 2 + a 2 (t)S l:j dx l dx 3 , (2) 
where a(t) is the scale factor. The Friedmann equation reads 

where the Hubble rate is H — a/a, dot denotes derivative with respect to cosmic time t and Mpi denotes the reduced 
Planck mass. We will assume that V((f>) is dominated by the vacuum energy Vo, therefore the previous equation 
simplifies to 



H 2 « (4) 
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The equation of motion for the scalar field reads 

4>{t) + 3Hj>(t) + m\m [1 + A6(fo - 0)] = 0, (5) 

where 9(cf>) denotes the Heaviside step function. 
The slow roll parameters are defined as 

With the assumption of a vacuum dominated inflationary model, the solution of Eq. ([5]) before to is 

<f> b (t) = <p ob U+(t) + 4> ob Ut(t), (7) 
where </>ob and (j>ob are the initial conditions, which we choose such that 4>ob is equal to zero. U^(t) are defined as 



Ut{t) = e\ H \ (8) 



with A^ defined as 



where /i 2 is /x 2 = If is small then A^ ps — ^- and A b s» — 3 (1 — ^g-J . After i the solution of Eq. (O is 

0a («) = 0OaC/ a + W + ^OaL/-(t), (10) 

where 0oa and </>oa are integration constants and the functions U^(t) are defined as 

ut(t) = ^i m , (ii) 

with Aj defined as 



where /Lt 2 is /^ 2 = /i 2 (l + A). U + (t) describes the slow-roll solution while U~{t) the rapidly decaying solution. 
Imposing the continuity of <p and <fi across to determines the integration constants <po a and <po a as 



i , A- -A 

t>0a — 90b 



L^o(A+-A+) ; ^ = 0()b ^L^ e ^o(A+-A-)^ 



(13) 



A a — Aj A a — \a 

In the plots of the following sections, we will make the following choice for the parameters of the model 

= 6 x KT 9 Af Pi , H = 2 x 10~ 7 M Pi , cf) 0b = 10M Pl . (14) 

We should note that the use of trans-Planckian field values might not be well justified in the context of the effective 
field theory, where higher dimensional operators are expected to destroy the flatness of the potential. 

In Fig. [TJ we plot we slow-roll parameters e, rj and for an amplitude of the transition of A = 2. We chose 
to = —1/Hln.H (which corresponds to conformal time r = —1). One can see that for our choice of parameters 
e is always small, while rj and -ft become large for some time after the transition time to- During this time the 
slow-roll expansion breaks down and this will introduce some scale-dependence in the power spectrum and produce 
non-Gaussian perturbations as we shall see in the next two sections. 
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FIG. 1: Plots of e (left), r\ (center) and jf^r (right) for A = 2. The heights of the peaks in the central and right plots are 
approximately 6A and —18.4(1 + A), respectively. While e remains always small, both rj and -^^j- become large after the 
transition and will source a large bispectrum. 



III. PERTURBATIONS 



In this section, we will discuss linear perturbations of the previous background. We shall present an analytical 
approximation to the mode function that will allows to obtain an analytical and accurate approximation to the power 
spectrum on scales sufficiently different from the scale fco set by the transition. Finally we integrate numerically the 
equation of motion of the curvature perturbation and obtain the power spectrum which we then compare with the 
analytical approximation. 

On the comoving time-slices, the scalar field fluctuations vanish, 8cj> — 0, and the three-dimensional spatial metric 
hij is perturbed as |l4j 

/^=a 2 e 27 %, (15) 

where tensor perturbations have been neglected because they do not contribute for the tree-level scalar bispectrum 
and 1Z denotes the curvature perturbation on comoving slices. 
The linear equation of motion for 1Z is 

4- (a 3 e^-n) -aed^ ^r- . 11 = 0. (16) 
dt V dt J dx l dx-> v ; 



In Fourier space and using the variable z = av2e the previous equation can be written as 

n'l + 2-n' k + k 2 n k = o, (17) 

z 

where prime denotes derivative with respect to conformal time t and k denotes the comoving wavenumber. 

Because we assume the potential is vacuum energy dominated, H is constant and the scale factor may be approxi- 
mated by that of a pure de Sitter universe a = —1/(Ht) and —Ht = ln(— Ht). We choose a(t — 0) = 1. The inflaton 
perturbation on flat hypersurfaces 6<j> is related to TZ at first order as 

54> = -^n. (is) 

The equation of motion for the field perturbation can be written as 

u" + (k 2 -—\u = 0, (19) 



where the variable u is defined as u = a5<fi = — sign(0)Afp;z7?.. At the time of the transition t = t$ (or t = to), Eq.([5]) 
implies that <p" is discontinuous which in turn implies that z" contains a Dirac delta function. Using Eq. p^|) one 
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can immediately see that also u" contains a Dirac delta function and u' is discontinuous at r = tq. One can evaluate 
this discontinuity as 

r-ra+i II fT +S illl (h(rn) 

D = lim / Z -dr = lim / t-dr = - m ^( To) |pi. (20 ) 

^ JTo-i Z ^°Jto-Z <P ( f ) ( T 0) 

Thus the matching conditions for u at the point r = To are given by 

^o + ) = "(To), u'(t+) = u'(t ~) + D «(r -), (21) 

where the superscripts — and + mean the lhs and rhs limits at the point To respectively. In terms of the comoving 
curvature perturbation TZ, these matching conditions imply the continuity of both TZ and TZ' at r = To. This is 
obviously consistent with the equation of motion for TZ, Eq. (jXTJ) , in which there is no Dirac delta function. In the 
next section, it is this last equation (fTT|) that we will integrate numerically to find the solution for the mode functions. 
Regarding the initial conditions we will assume the standard Bunch-Davies vacuum for 8(j) at r — > — oo. At sufficiently 
early times, r — > — oo, the mode functions are well approximated by 

e -ifcr / ; \ 

ut = vs vmV-Vr)' < 22) 

where we denote the mode function at r < To by Uf,. In our numerical analysis, we will use this previous equation to 
set the initial conditions at the initial integration time Tj for each mode u when it is inside the horizon and for r < To. 
Then we change variables to TZ and integrate numerically Eq. (|17[) . 

For long wavelength modes compared with fco, i.e. fc < fco, the transition in the mass happens when the modes are 
outside the horizon and Eq. (|22[) is a good approximation until the modes cross the horizon. So we will use it to set 
the initial conditions at some time before horizon crossing and then solve Eq. (|16[) numerically. On the other hand, 
for modes with wavenumber greater than fco, the sudden change in the mass happens before horizon crossing and we 
will set the initial conditions using Eq. (|22[) at the time t = Ti < To when it is still a good approximation for the 
mode function. 

To quantize the curvature perturbation, we follow the standard procedure in quantum field theory. We promote TZ 
to an operator that is expanded in terms of creation and annihilation operators and mode functions as 

TZ(t, k) = TZ(t, k)a(k) + TZ*(t, -k)a f (-k), (23) 

where a and a) satisfy the usual commutation relation [a(k), a^(k')] = (27r) 3 (^ 3 )(k — k'). 
The power spectrum of the curvature perturbation is given by 

2r^ /2 (fc) =V2tf\Tl k {tf)\, (24) 

where tf is the time at which inflation ends. Piz(k) has been measured observationally to be -pR.(fco) = 2.95 x 
l(T 9 yl(fc ), where A(k ) is A(k ) ~ 0.71 - 0.75 and fc = 0.002Mpc _1 [HI]. 



A. Analytical approximation to the mode function 



In this subsection, we will present an analytic al ap proximation to the mode function. This approximation has been 
discussed in previous work, see for example (97l . IllOl ] . 

For modes with k < fco, the transition happens when the modes are outside the horizon and the only possible time 
variation of TZ is through the discontinuity of z' / z present in Eq. (|17p. however z' j z = az/z which means that it 
grow s exponentially with time and so the effect of the discontinuity become more and more irrelevant. In fact, Ref. 
[1 191 ) has shown that the only expected modification with respect to the usual result is for scales around fc = fco- 
Therefore, for fc < fco an d for r < t^, where t& is the horizon crossing time for that mode, we will approximate the 
mode function by Eq. (|22|) . 

For modes with wavenumber fc > fco, Eq. (|22j) is still a good approximation for t < to and we shall use this fact. 
For some time after the transition, the slow-roll expansion breaks down and one can not use the slow-rolling solution 
for the background, however for sufficiently large fc, the mode is well inside the horizon at the time of the transition 
and this small violation of slow-roll becomes negligible because the mode equation (|T9)) is dominated by the term 
proportional to fc 2 and the slow-roll violating term z" / z can be safely neglected. Hence, Eq. (|2"2"j) is still a good 
approximation even for t > t , including the short non-slow-roll period after the transition. However, the solution 
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for the mode function u a after r = tq is no longer given by just the positive frequency mode function v, instead it 
will be a mixture of positive and negative frequency modes as 

u a = a k v + (3 k v*, (25) 

where the Bogoliubov coefficients a k and (3 k can be expressed in terms of v, v* and the outgoing mode function u a as 

a k = -i (v*'u a - v*u' a ) , (5 k = i (v'u a - vu' a ) , (26) 

and they satisfy \a k \ 2 — \/3 k \ 2 = 1. The matching conditions (|21|) imply 

u a{To) = v , <(r^) = v' Q + D Q v , (27) 

where vq = v(t ) and together with Eqs. (|26|) applied at the time t = they can be used to find the Bogoliubov 
coefficients as 



a k = 1 + iD v Vq = 1 + i— - ( 1 + j- 1 , 2 

2fc V (fcror 



/3 fc = 



• n „ 2 ■ ^0 



: -Do 
2fc 



fcr 



— 2ikro 



(28) 



In terms of the original variable, the comoving curvature perturbation, the analytical approximation for the mode 
function that we use is 



TZ(t, fc) 



Mpia{T) 



v(r.k) 
v(r.k) 

a(k)v(T,k)+P(k)v'(T,k) 
a(k)v{T,k)+l3(k)v' (r,k) 



k < kg and r < r k 
k < ko and r > r k 
k > ko and t < tq 
k > kg and tq < t < r k 
k > kg and r > T k 



(29) 



where the difference between the last two branches and between the second branch and the first and the third 
branches is that we fixed the value of the slow-roll parameter e to its horizon crossing value, i.e. when the modes 
are super-horizon. We have assumed that sign(0) = — 1. Note that Tt(r,k) is discontinuous on the line k = ko 
and for t > tq. The time derivative is discontinuous on the horizon crossing line kr k = —1, in particular the time 
derivative is discontinuous at the point (to, ko), i.e. 7^' (to, fcg~) ^ 7?.'(r , k^), we choose H'(tq, k ) = H'(t ~ , Note 
that TZ'(tq , ky) = TZ'(to, k^). We have solved numerically the mode function equation of motion (fT6|) an used the 
analytical solutions for the background given in section |n] and we found that the previous analytical approximation 
(1291) is a good approximation for the mode functions at any time r and for scales far from ko even if the amplitude of 
the transition A is of order one. As expected this approximation deteriorates for scales near fco, however for A small 
this continues to be a good approximation. 



B. The power spectrum of the primordial curvature perturbation 



In this subsection, we shall present the result for the power spectrum calculated numerically and compare it with 
the power spectrum given by the analytical approximation. 

As discussed previously, for modes with wavenumber k such that k -C fco there can be no super-horizon evolution 
of 1Z and the curvature perturbation power spectrum has to be given by the standard formula, which in our present 



27r|0(t fc )| 2TA+0ofc V * ' 



where t k is the horizon crossing time, fc = a(t k )H. This is the same equation as in the case of no transition. Our 
analytical approximation consists of using this result all the way down to the scale fco however for scales near fco this 
might not be a good approximation. 

For modes with fc > fc , using Eq. (|2"5)) in the power spectrum definition (fM| one can find that the spectrum is 

P ^ /2(fc) = ^^ |afc - /3fe| ' (31) 
2it\4>{tk)\ 
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FIG. 2: Plots of P TC (fc) for A — 0.01 (left) and A = 2 (right). The red dashed line is the numerical result and the green 
continuous line is the analytical approximation. It can be seen that the analytical approximation is good for scales different 
from fco and in particular it provides a good approximation for the small-scale damped oscillations. On the lhs plot and for 
fc > lOfco the numerical result deviates from the analytical approximation by at most O(0.5%), this is due to the numerical 
error of our integrator. For large A the amplitude of the power spectrum on small scales is significantly different from the 
amplitude on large scales. 



where the Bogoliubov coefficients are given in Eq. (|28|) . In the following plots of our analytical approximation we will 
use the previous result for scales k > k but one should keep in mind that in principle this is expected to be a good 
approximation only for scales k 3> fco- 

For fc ^ fco, one can find a simple expression for the power spectrum |lldj | 



H 2 



27r|0(t fe )| 



For later use we define P> /2 (fc) as 



Pi /2 (fc) = 



H 



27rAj0o6 



§(l + cos(2fc T0)+O (-i- 



P< /2 (k) 
1 + A ' 



1/2 



(32) 



(33) 



This is the asymptotic value of the power spectrum for scales much smaller than the transition scale fco- We should 
note that if A is large then the amplitude of the power spectrum for scales much smaller than fco is significantly 
different from the amplitude for scales much larger than fco. The ratio of these amplitudes is proportional to f + A. 
The spectral index is defined as 



1 



dlnPn(k) 
din k 



= -2e 



V ; 



0( M 4 ) > 0, 



(34) 



where in the last equality we particularized the result for our model and we used the formula for Pji given on large 
scales Eq. (|30p . This means that for scales larger than fco the spectral index in the present model is positive. For 
scales smaller than fco the power spectrum contains oscillations which are damped. For fc ^> fco, the spectral index 
approaches a constant given by 2/3jz?. Recent observations favor a negative spectral index [7j, however a scale 
invariant spectrum is still allowed [l20(, this constraints the parameters /i 2 and /i 2 of the model to be much smaller 
than one. 

As can be seen from Fig. [2] the analytical approximation agrees well with the numerical results at large and small 
scales compared with fc . As expected, this approximation loses some accuracy for modes around fco- This is because 
these modes leave the horizon around the time of the transition, for them the effect of the violation of the slow-roll 
conditions is non-negligible and the approximate solution (|22p is no longer good. 
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IV. THE BISPECTRUM 



In this section, we will present th e calculation of the bispectrum of the primordial curvature perturbation using the 
so-called in-in formalism jl2ll Il22| | . We will use the analytical approximation of the mode function described in the 
previous section and therefore we will be able to obtain analytical expressions for the leading order bispectrum for 
the present model in certain interesting limits. 

In order to use the machinery of the in-in formalism to compute the tree-level three-point correlation function (or 
bispectrum) one needs to calculate the cubic-order interaction Hamiltonian, see for example [76] for a review about 
this procedure. 

The third order action after ignoring many total derivative terms that appear when one simplifies the acti on by 
using integrations by parts has been know since the seminal work by Maldacena [HJ and can be also found in [26L 1 123l j . 
it reads 



Sa = M 



pi 



dtd 3 x 



a 3 e 2 KK 2 



a 3 e dr\ r 



2 K{dK) 2 - 2ae1Z(d1Z)(dx) 



~—^TZ 2 TZ + (dK)(d X )d 2 x + -(d 2 K)(d X ) 2 + 2 
2 at 2a 4a 



\ 5L 




/ sn 


1- 



where we should note that no slow-roll approximation has been made and we define 

fdd 2 X 

= a 

l 



X 



a 2 ed- 2 K, 



SL 
STZ 



Hd 2 X - td 2 K 



1 



f(K) = -K1Z 



1 



4a 2 H 2 



V dt 

\-(dlZ){d1Z) + d^ididjidiTZdjU))] 



(35) 



(36) 



1 



2a 2 H 



[(dTZ)(d X ) ~ d-^didjidiTZdjx))] ■ (37) 



In the previous equations, d~ 2 denotes the inverse Laplacian and 6L/SlZ\i denotes the variation of the quadratic 
action with respect to the perturbation 1Z. In order to proceed, usually one performs a fi eld redefinition to eliminate 
the terms proportional to the linear equations of motion. Recently, it was shown [124] that instead of doing field 
redefinitions one obtains the same result by working with the original variable 1Z and without ignoring the boundary 
terms that were omitted in Eq. (|35l) . In the present work we will follow a mixed procedure, i.e. we will use a field 
redefinition as 

K^K n + f(K n ) (38) 
to el iminate the last term of the third order action (|35|) and we will keep the boundary terms as in [124] (see also 



[123). 

The previous field redefinition involves terms that contain at least one derivative on 1Z therefore if one chooses to 
evaluate them after horizon crossing they should give a negligible contribution to the bispectrum. This implies that 
at such time the bispectra of 7Z and lZ n are equal. For this reason from now on we will drop the subscript n and 
identify lZ n with 1Z. 

In Maldacena's calculation [l4| and in the following ones [26|, Il23| the first three terms in the second line of Eq. 
(|35[) which are higher-order in the slow-roll expansion were properly neglected because these authors work at leading 
order in slow-roll. In the present case, e is small but r\ and rf may be large, this means that the dominant contribution 
to the three-point function comes from the following terms in the total third order action: 



S 3 D M 2 Pl 



dtd 3 i 



a 3 e dr\ 



SL 



— — ^TClZ + 2±1Z Z — 



2 dt 



51Z 



M 



pi 



dtd i x— 
dt 



- ^n 2 d 2 x 

2 



(39) 



where the last total time derivative term comes from the boundary interaction terms recently obtained in [126| (see 
also 124]). Instead of using the previous third order action we found it is more convenient to do one integration by 
parts in time to simplify the action to 



S 3 D M 2 Pl I dtd 3 x 



er]a 3 1Z1Z 



32 



er) 



a1Z 2 d 2 1Z 



(40) 



This form of the action makes it clear that TZ is constant outside the horizon to non-linear order as expected . 



1 This procedure is equivalent to the usual one used in for instance Refs. l67j |68l| . 
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FIG. 3: Plots of e (left) and 77 (right) for A — 2. The blue dashed lines are the plots of the approximated expressions Eqs. (|44[) . 
The red continuous lines are the plots of the exact analytical solutions found in Sec. [TT| The height of the peak in the rhs plot 
is approximately 6A. While e remains always small, 77 becomes large after the transition and will source a large bispectrum. 
The simple analytical approximation to r\ becomes poor some time after the transition but from this time r\ is small anyway 
and we expect a very small error to the bispectrum coming from this deviation. 



The interaction Hamiltonian in conformal time is 



d 3 s 



e^aRK' 2 + ^-aTZ 2 d 2 TZ 



H int (r) = M 2 pl J 

The tree-level three-point correlation function at the time r e after horizon exit is 

(0|^(re,ki)^(Te,k 2 )^(T e ,k 3 )|0) = -* / dra(0\ [Ufa, ki)&(r e , k 2 )7t(r e) k 3 ), H int {T 



(41) 



(42) 



where |f2) and |0) denotes the interacting vacuum and the free theory vacuum respectively. [,] denotes the standard 
commutator and the interaction Hamiltonian Hj n t is used to evolve the free theory vacuum to the interaction vacuum 
at the time the three-point function is evaluated 14]. 
More explicitly, the bispectrum is 



(fi|&(0,k 1 )&(0,ka)£(0,k3)|fi> « (27r) 3 ,5 (3) (ki + k 2 + k 3 )2M 



2 o 



pi 



7e(0,ki)7e(0,k 2 )7e(0,k 3 ) / dTr)ea^TZ*(T, ki) x 



{2W* (r, k 2 )W* (r, kg) - k 2 TV (r, k 2 )TZ* (r, k 3 )) 



two perms., 



(43) 



where we have set r e = and we evaluate the integral from To instead of — 00 because for r < To, slow-roll is valid, 
the slow-roll parameters are small and therefore the integral is small. We expect that the approximation of setting 
r e = in the integral is good enough because some time after horizon crossing the mode function becomes constant 
and the integral approximately vanishes, "two perms." means the other two permutations of ki,k 2 and k 3 . 

In the following calculation of the bispectrum, we will approximate the exact analytical solutions for the slow-roll 
parameters found in Sec. Ullas 



e(r) 



4>l b {K? 

2M 2 



a{rf 



L Pl 



1 - 



.4 



1 + A Vto 



77(7 



-6 



l-l±A(m r 



(44) 



which are a goo d approximation to the exact analytical expressions for times after to, as can be seen in Fig. [3] 

Following [671 ] . the meaningful quantity to plot to show the shape of the bispectrum when the three wavenumbers 
k^ i = 1,2,3 have comparable magnitudes is 



g(fci,fc 2 ,fc 3 ;fc») 
kik 2 k 3 



(hk 2 k 3 ) 2 



<5( 3 ) (kr+k 2 +k 3 ) (2Tr) 7 p2(fc„) 



(Q|ft(T e ,k0ft(T e ,k 2 )ft(T e ,k 3 )|n), 



(45) 
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where fc* denotes a pivot scale at which the power spectrum P-jz is to be evaluated. In the previous equation one 
should use the formula for the power spectrum amplitude as given by Eq. ([50")) or by Eq. ([3"5"f depending on whether 
one is interested in large or small scales respectively. In this way the oscillations seen in the following plots are truly 
bispectrum oscillations and are not due to the oscillations already present in the power spectrum. If this quantity 
is at least of order of a few then there is hope that for instance the Planck satellite will measure it. In the present 
work, because the asymptotic amplitudes of the power spectrum for scales much larger and much smaller than the 
transition scale can be very different if A is large, it makes sense to define two separate functions Q < and Q > to study 
the bispectrum when the three wavenumbers have magnitudes smaller and larger than fco respectively. The two new 
functions are defined as 



G<(k ll k 2 ,k 3 ;h) _ 1 (kxk 2 k 3 ) 



>2 



k x k 2 k 3 <J(3) (ki + k 2 + k 3 ) (2tt) 7 P| [h 



■<0|^(T e ,k 1 )^(Te,k 2 )^(T e ,k 3 )|0), (46) 



G>(ki,k2,ka;k*) 1 (fcifc 2 fc. 



?>) 



\2 



■(0|^(re,k 1 )^(re,k 2 )^(r e) k 3 )|n}. (47) 



k X k 2 k 3 ,5(3) (k! + k 2 + k 3 ) {2vfPl{K 

For comparison the so-called local model of the bispectrum is 

(^(k!)^(k 2 )^(k 3 )) iocai = (2^) 7 <5 (3) (ki + k 2 + k 3 ) A /jViP 2 (4g) 

where the sign of ,f/y/, was chosen to agree with the WMAP team definition Q but it has the opposite sign of the 
definition present in [14 [67 1 for example. We should stress that the definition of Jnl in the present paper is different 



from the definition in [7| but both definitions agree for a scale invariant power spectrum. 

Following the definition of the function Fnl, which is useful to study the squeezed limit of the bispectrum and was 
first introduced in Ref. [70|, here and for the reason discussed above, we also define two functions F^ L and Ff^ L as 

, , , lOfci k 2 k 3 (fci , k 2 , k 3 \ fc*) 

F< L (k Uk2 ,k M , 3 £y <l ;; fc y (49) 



t ; . , _l0k 1 k 2 k 3 G > (k 1 ,k 2 ,k 3 ;h) 



F NL {k 1 MM;K) = ^ T - w ^ , (50) 

which reduce respectively to G<(ki, k 2 ,k 3 ; fc*) / '(k\k 2 k 3 ) and G>{k±, k 2 , k 3 ; fc*)/(feife 2 fc 3 ) in the approximately equilat- 
eral case and to Jnl for the local model. 



A. The squeezed limit 



Using the analytical approximations for the mode function ([2TR) and for the slow-roll parameters (|4"4")l one can 
calculate the integrals analytically, for scales larger than ka, or numerically for any scale. 
In the squeezed limit for large scales, i.e. k\ <C k 2 , k 3 ~ k and k <^ ko, 



1 ^(8 + 3^1) 
12 (1+.4) 2 



A+ 



(r )a 2 (T*) 



2 A 



a(T kl )a 2 {T k ) 



(51) 



The suppression of -F/vx by the square of the ratio fc/fco has b een previously found using the next-order gradient 
expansion method [73| and using the consistency relation |l27| . In this limit and for values of A of order one the 
quantity Fnl is too small as can be seen in Fig. 01 

To calculate the bispectrum on small scales we used a further approximation in the analytical expressions for the 
mode function Eq. (|2T)1) . We fixed the value of e(r) to the horizon crossing value of that mode e(Tfc) for any time after 



11 



10 20 30 40 50 



NL -0.3- 



FIG. 4: Plots of F< L (fc /500, ko/n, ko/n; k ) for A = 2. For the red line, the bispectrum integral in Eq. (|43[) was calculated 
numerically and we used the analytical approximation for the mode function. The green and black lines are the plots made 
using the first and second parts of Eq. (|51f> respectively. It can be seen that the lines are practically indistinguishable. For 
scales k ~ ko, i.e. k of order one, our analytical approximation for the mode functions should break down and the result in 
this plot should not be taken seriously. 



and before horizon crossing. We found for the squeezed limit on small scales, i.e. k% <C fc 2 , fc 3 ~ k and ko <C fci 



Fj\n,{k\,k,k;k^) ~ I — r 



XaJ a 6A ^(r ) V W ( T k) 
2k + k 1 k 1 . f 2k + fci 



2,V 



5A 



A 



h k \ k Q 
2k + fci fci . / 2fc + fci 



, 2fc + fci \ „ , f-2k + ki 
+ (2 -3fA) cos — -— — -3/Acos' 



4(1 + A) 4 



ko 



■ sm 



fco 



ko J 
f 2k + fci 

(2 - 3/A) cos f — - — - ) :-!./•. I .s 



-2A; + fci 
fco 



(52) 



where the parameter / should be equal to one. If / = this is equivalent to set the Bogoliubov coefficients to 
= 1,/3/c = (the expected behavior in the infinite frequency limit). We found that the previous additional 
approximation of fixing the value of e(r) to the horizon crossing value of that mode e(rfe) for any time after and before 
horizon crossing is not a good approximation in the calculation of the amplitude of the bispectrum, this is because 
for some time after the transition the value of e(r) varies significantly due to the presence of the background decaying 
mode of the scalar field. However, as can be seen in Figs. [5] and [6] this approximation reproduces the shape very 
well. If one allows e(r) to vary, the amplitude turns out to be larger than in the fixed e(r) case. In the varying e(r) 
case, we were unable to calculate the bispectrum integral analytically so we proceeded numerically in order to find 
the amplitude, the result for A — 2 can be found in Figs. Elands We were able to find a simple and accurate fitting 
formula to the numerical result for the amplitude of F^ L on small scales as 



F^ L (50fc , lOOOfco, 1000fco; fco) 



20(14(1 + A) 



CiA c »e c > A , 



(53) 



where the constants are 



Ci « 2873.31, C 2 2.00, C 3 ~ 0.02 or C x « 2961.32, C 2 = 2, C 3 = 0. 



(54) 



In the range of values of A as A — [0.01, 10], for the first and second set of constants the approximation is good to 
2% and 15% respectively. In Fig. [7] we plot the value of F^ L (50fco, lOOOfco, lOOOfco; fco) as a function of A using the 



numerical result and the above fitting formula (|53|) . 

As can be seen in Eq. (|52| or in Fig. [5j the envelope of the sinusoidal grows with fci, this can have important 
consequences regarding the magnitude of the non-Gaussianity of this model and its potential observability even for 



12 




FIG. 5: Plots of F^ L (nko, lOOOfco, lOOOfco; fco) for A = 2. For the red dots, the bispectrum integral in Eq. (I43|) was calculated 
numerically and we used the analytical approximation for the mode function. The black line is the plot made using the last 
line of Eq. (|52[) . The blue dashed line is the plot made using only the first term in the last line of Eq. (|52[1 . It can be seen 
that both lines are almost indistinguishable and are a good approximation for the red points. Both black and blue lines were 
normalized at the point k = 50 with the normalization fitting formula Eq. (|53|) and using the first set of constants of Eq. (I54[l . 




FIG. 6: Plots of F^ L (50ko, nko, nko; ko) for A — 2. For the red dots, the bispectrum integral in Eq. ([43j) was calculated 
numerically and we used the analytical approximation for the mode function. The black line is the plot made using the last 
line of Eq. (|52|l . The blue dashed line is the plot made using only the first term in the last line of Eq. (|52[) . The two lines are 
almost indistinguishable in the plot and they are good fits to the red points. Both black and blue lines were normalized at the 
point k — 1000 with the normalization fitting formula Eq. (|53l) and using the first set of constants of Eq. (|54ll . 

A relatively small. This growth might be worrisome because the bispectrum would increase without bound for very 
small scales, however one should note that at some scale our toy model, i.e. replacing an otherwise smooth transition 
by a Heaviside function, ceases to ma ke sense. This scale gives a natural cut-off below which the previous result can 
not be applied. For instance in Refs. jllOl Il28j it was argued that for the particular model under consideration there 
the Heaviside function approximation was valid for k < 10 3 fco- For scales smaller than this cut-off and if the width of 
the transition is taken into account then the linear growth should stop and the amplitude should decay to small values 
because for these sufficiently high frequencies the transition becomes an adiabatic change of the slow-roll parameters 
and in this case we expect nearly Gaussian perturbations. 

With a simple order of magnitude argument it is possible to estimate the range of wavenumbers Ak over which 
we expect to have a large deviation from Gaussianity as follows. If instead of a sharp transition we have a smooth 
transition with a width in field space given by A0 then for the vacuum dominated potential Eq. ([1]) one can 
find that the transition happens over a number of e-foldings AiV as AiV ~ (3A0)/(/i 2 </>o). Using the fact that 
Afc/fco ~ To/ At ~ I /(AN) one estimates that the range of scales that can be affected by the transition is Ak ~ 
ko/(AN) ~ (fi 2 (f> )/ (3A(f))ko. One can easily see that for a sharp Heaviside transition A(f> = then Afe — > oo. 
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FIG. 7: Plots of F^ L (50ko, lOOOfco, lOOOfco; ko) as a function of A. For the red dots, the bispectrum integral in Eq. (|43[) was 
calculated numerically and we used the analytical approximation for the mode function. The continuous green line and the 
black dashed line are the plots made using the fitting formula Eq. ()53[) with the first and second set of constants as in Eq. (|54[) 
respectively. It can be seen that the green line fits well the red points. The line deviates from the red points at most 2%. 



— — -0.6 



FIG. 8: Plots of Q< (k, k, k; ko)/k 3 with k — ka/n for A — 2. For the red line, the bispectrum integral in Eq. (|43p was calculated 
numerically and we used the analytical approximation for the mode function. The green and black lines are the plots made 
using the first and second parts of Eq. ((55} respectively. It can be seen that the lines are practically indistinguishable. For 
scales k ~ ko, i.e. k of order one, our analytical approximation for the mode functions should break down and again the result 
in this plot should not be taken seriously. 



B. The equilateral limit 

In the equilateral limit and on large scales, i.e. if k <C ko, we find 



G<(k, k, fc; k*) 
P 



27 A(8 + 3A) /A+ 
160 (1 + A) 2 



A + 

2 / 7 \ 2 / T¥ , \ o / 



a x >> (T )a 2 (T* ) 
a 3 (T fe ) 



2A^ 



(55) 



The plot of Q < (k,k,k;ko)/k 3 for A = 2 can be seen in Fig. [8] The amplitude of the graph is too small to be 
interesting observationally. This plot highlights the fact that it is very difficult to generate large non-Gaussianity in 
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single-field models of inflation when the scales of interest are well outside the horizon. This is a direct consequence of 
the constancy of the curvature perturbation on these scales and of the fact that at horizon crossing the perturbation 
was very Gaussian because the transition had not happened yet. 

In the equilateral limit and on small scales, i.e. if fco <C k, and using the extra approximation in the analytical 
expression for the mode function Eq. (|29|) of fixing the value of e(r) to the horizon crossing value for any time, we 
find 



G>(k, fc, fc; fc* 
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,4A 



" ( T *) ( a 4 (r ) 
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cos 



9Af ( fc 
— — cos — 
2 \k 



(56) 



where the parameter / should be equal to one. If / = this is equivalent to set the Bogoliubov coefficients to 
afc = l,Pk = (the expected behavior in the infinite frequency limit). The term with the linear growth in the 
previous equation, i.e. the term proportional to k/ko, remains even if / = 0. This means its origin is not from 
the mixture with the negative frequency modes. In fact, this term is present even in the standard slow-roll scenario 
if we artificially switch off the interaction for t < To, switch it on at tq and perform the integration from tq. As 
argued before, if the finite width of a realistic transition is taken into account then for scales smaller than the scale 
set by this width we expect the previous result to be unapplicable and the amplitude of Q> /fc 3 should decay to small 
unobservable values. 

The plot of the previous expression for A — 2 is displayed in Fig. |H1 As in the squeezed limit, the extra approximation 
of fixing the value of e(r) to the horizon crossing value is not a good approximation in the calculation of the amplitude 
of the bispectrum as given by the previous equation. However, to obtain the correct shape that approximation is 
enough as can be seen in Fig. [5] In order to obtain the correct amplitude of the bispectrum we calculated the integral 
numerically, but using the analytical approximations for the mode functions, and we were able to find a simple and 
accurate fitting formula to the numerical result for the amplitude of G> / fc 3 on small scales as 



G>(k, fc, fc; fc ) 



fc :! 
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k=100k 



200,4(1 + A) 



CxA e *e G * A , 



(57) 



where the constants are 



C x 



-5157.87, C 2 2.01, C 3 



0.03 



C x 



-5399.03, C 2 = 2, C 3 = 0. 



(58) 



In the range of values of A as A = [0.01, 15], for the first and second set of constants the approximation is good to 
5% and 50% respectively. In Fig. [TU]we plot the value of G>/k 3 with fc = 100fco and fc* = fc as a function of A using 
the numerical result and the above fitting formula (I57|) . 

It is worth mentioning that in Ref. [68j, the authors proposed an approximate sharp feature ansatz to be used in 
data analysis as 



fNL{k) sin(K/k + phase), 



(59) 



where K = fci + fc 2 + k% and ^Nh{k) is some envelope. This form seems to fit well with the current model. Finally, it 
is important to mention that the linear growth with k that we found in the present model (an Heaviside step function 
in the inflaton's mass) is in principle observationally distinguishable from the model considered in jl03l | where it is 
the slow-roll parameter e that has a Heaviside step function discontinuity and the author found a quadratic growth 
in K 2 . 



V. CONCLUSION 



In this paper, we introduced an inflationary model driven by a single standard kinetic term scalar field. The field's 
potential is dominated by the vacuum energy and the mass term changes abruptly at a point. This abrupt change 
is modeled by an Heaviside step function. This might be seen as a simple toy model for a first order mass phase 
transition occurring during inflation. 



2 We would like to thank Xingang Chen for pointing these out to us. 
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FIG. 9: Plots of S> (k, fc, fc; fco)/fe 3 with k — kUq for A = 2. For the red dots, the bispectrum integral in Eq. (|43p was calculated 
numerically and we used the analytical approximation for the mode function. The black line is the plot made using the last 
line of Eq. (|56|) . It can be seen that the black line is a good approximation for the red points. The black line was normalized 
at the point k = 100 with the normalization fitting formula Eq. (|57l) and using the first set of constants of Eq. (|58[) . 

A 

0.01 0.05 0.1 0.5 1 5 10 



J I I L 




FIG. 10: Plots of g > {k,k 1 k;k )/k 3 with k = lOOfco as a function of A. For the red dots, the bispectrum integral in Eq. f!43[) 
was calculated numerically and we used the analytical approximation for the mode function. The continuous green line and 
the black dashed line are the plots made using the fitting formula Eq. (|57[l with the first and second set of constants as in Eq. 
(|58[1 respectively. It can be seen that the green line fits well the red points. The line deviates from the red points at most 5%. 

Under the vacuum domination assumption, we have solved the background evolution analytically. As expected for 
large step sizes and for some time after the transition the slow-roll approximation breaks down. 

We have considered linear perturbations and presented an analytical approximation for the mode function that for 
scales different from the step scale kg is a good approximation to the numerical mode function for any time. We have 
computed the power spectrum of the curvature perturbation and showed that for scales smaller than ko it contains 
damped oscillations with "angular frequency" 2/fco- This behavior is captured well by our analytical approximation 
of the spectrum. 

Models with temporary violations of slow-roll are expected to produce large non-Gaussianity. For our choice of 
parameters, we showed that the dominant contribution to the bispectrum comes from the terms containing 77 in the 
third-order interaction Hamiltonian. Using the analytical approximation of the mode function, we have computed 
(analytically when it was possible, otherwise numerically) the bispectrum produced by these terms in certain inter- 
esting limits. We analyzed two well known limiting triangle configurations of the three momentum vectors on which 
the bispectrum depends on. 
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First, for the equilateral limit, we found that for scales much larger than ko the quantity G</k 3 decreases with the 
square of the ratio k/ko which implies that for values of A of order one, its value is too small to be of observational 
interest. For scales much smaller than ko it oscillates with an "angular frequency" of approximately 3/ko- On these 
scales Q y /k 3 is linearly enhanced towards large k as G>/k 3 oc k/ko- This enhancement may be used to push the 
value of G>/k 3 to within observational range even for a small amplitude A of the step. It would be interesting to 
consider observational constraints on this kind of strongly scale dependent and oscillatory bispectrum models but this 
is outside the scope of the present paper. 

Secondly, we considered the squeezed limit of the bispectrum. For scales much smaller than ko, again we found an 
oscillatory behavior in k with "angular frequency" 2/fco- When the size of the smallest side of the squeezed triangle 
is k\ 3> ko then the quantity F^ L has as enhancement factor of k\/ko- This might have important consequences 
regarding the detectability of the present signal in the data. For larger scales, we found that the quantity F^ L is 
suppressed by the square of the ratio k/ko which implies that for values of A of order one, its value is too small to be 
observational relevant. 

We have shown that for a wide choice of the model parameter A, the squeezed and equilateral limits of the bispectrum 
can be large for scales much smaller than ko. This represents a significant enhancement with respect to the canonical 
single field slow-roll model result. In particular, the squeezed limit of the bispectrum can be large in comparison 
with any single field model of inflation known to the authors. See however Refs. [741 l75l for non- vacuum initial state 
scenarios which also give significant enhancement factors but which have been shown 75[ not to be enough to put 
the CMB bispectrum within the reach of the Planck satellite. It would be very interesting to repeat that analysis for 
the present model. 

The analytical approximation for the mode function used in this paper breaks down for scales close to the scale set 
by the transition so the present method cannot be applied to compute the bispectrum at these scales. In order to 
perform that computation one has to numerically integrate the equation of motion for the curvature perturbation and 
then calculate the bispectrum integral numerically. We leave this for future work. Finally, it would also be interesting 
to study the trispectrum of this model and this will be presented in a forthcoming publication. 
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